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Two zone equilibrium and kinetic free model proposed by the authors in their earlier work [Ratnadhariya 
JK, Channiwala SA. Two zone equilibrium and kinetic free modeling of gasifier. Proceedings of the 12th 
European Conference and Technical Exhibition on Biomass for Energy, Industry and Climate Protection. 
Amsterdam, The Netherlands; 2002. p. 813-816], offers gas composition, temperature profile and gasifier 
performance parameters for two zones. This model does not predict composition and temperature 
profile of pyrolysis zone, which is stated to be a precursor for gasification. Looking to this fact a three 
zone equilibrium and kinetic free model of biomass gasifier is proposed in the present work. In this three 
zone: first zone of the model is drying and pyrolysis zone combined together; second zone is oxidation 
zone; and the third zone is the reduction zone. Each zone has been formulated with: (i) reaction stoi¬ 
chiometry; (ii) constituent balance; and (iii) energy balance along with a few justifying assumptions. This 
model clearly provides an operating range of equivalence ratio and moisture content for the woody 
biomass materials. Further, this model facilitates the prediction of the maximum temperature in the 
oxidation zone of gasifier, which provides useful information for the design of the gasifier and selection 
of the material for the construction. The merits of the model lies in the fact that it is capable of handling 
predictions for all category of biomass materials with a wide operating range of equivalence ratio and 
moisture content in all of the three principal zones of the gasifier. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

A novel approach has recently been introduced by the authors to 
formulate a two zone equilibrium and kinetic free model of the 
gasifier based on a few justifying assumptions, which indirectly 
considers reaction kinetics through exothermicity and endo- 
thermicity of principal oxidation and reduction reactions [1]. The 
present model is an extension of the same effort. 

The present model identifies the sequence of principal reactions 
and a degree of each reaction is considered to be inversely 
proportional to their exothermicity or endothermicity and hence 
does not need any assumptions about thermodynamic equilibrium. 
It predicts composition and temperature levels in three zones 
simply from mass and energy balances. 

2. Physical description of model 

Fig. 1 gives the overall view of the physical model of the 
downdraft gasifier along with possible reactions in the different 
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zones. Drying and pyrolysis zones are clubbed together. The prod¬ 
ucts of pyrolysis zone, i.e., 1st zone are assumed to be char, CO, CO 2 , 
H 2 , H 2 O, CH 4 and C 2 H 2 [2-8]. Higher hydrocarbons are not 
considered as they crack during secondary pyrolysis in the gasifier 
as observed by Antal et al. and others [5,9-12]. Similarly, many 
recent studies on pyrolysis gives gas composition at different 
temperatures and different heating rates and their results clearly 
indicate release of CO, CO 2 , H 2 , H 2 O, CH 4 and tar with a higher 
concentration of lighter components [5-8,13]. 

The initial association of C, H and O with each other during 
formation of biomass [14-19] along with process parameters, i.e., 
temperature, rate of heating, etc., usually govern the product 
quality and yield in the pyrolysis zone [9-11,13,20,21]. Channiwala 
[14] has observed that the majority of HHV correlations are based 
upon the different assumptions regarding association of oxygen 
with hydrogen or carbon. Dulong [15] and Strache and Lant [16] 
have used the concept of available hydrogen and assumed fuel 
oxygen to be totally associated with hydrogen. Steuer [17] and 
Sumegi [18] have considered the association of fuel oxygen with 
both hydrogen and carbon while Vondrecek [19] has assumed more 
availability of hydrogen. Sumegi [18] has assumed (3/8) O to be 
associated with C and (1/16) O to be associated with the hydrogen 
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their reactions [2,3,27-29]. Finally, CO and CO 2 from the pyrolysis 
zone will be simply added up to the CO and CO 2 produced by char 
oxidation. CH 4 and C 2 H 2 are simply assumed to be carried forward 
from the pyrolysis zone to oxidation zone due to their lower 
burning velocities and oxygen deficiency [2,27-29]. Thus, the 
products of the oxidation zone are char, CO, C0 2 , H 2 0, CH 4 , C 2 H 2 
and N 2 [2-4]. The overall composition and temperature in this zone 
is obtained through mass and energy balance. 

The products from the oxidation zone enter reduction zone. 
The possible products of the reduction zone are assumed to be CO, 
CO 2 , H 2 , H 2 O, CH 4 , C 2 H 2 and N 2 [2,5-8]. The water-gas shift 
reaction is assumed to have a degree of reaction as a unity at the 
entry of reduction zone as it is a gas phase homogeneous reaction 
[2-4,32-37]. Char carry over is assumed to be zero [38,39]. CH 4 , 
C 2 H 2 and N 2 from the oxidation zone are carried forward to 
product gas [2,3,27-29,40]. Balanced carbon will be utilized by 
principal reduction reactions in inverse ratio of their endo- 
thermicity, i.e., Boudouard reaction will consume 43.22% carbon 
while water-gas reaction will consume 56.78% carbon [3,27-29]. 
The composition and temperature in this zone is obtained through 
mass and energy balance. 

3. Model formulation and solution 


| CHAR + ASH | CO, C0 2 , H 2 0, CH 4 , C 2 H 2 , N 2 


ZONE-III 

REDUCTION ZONE 

C + C0 2 —2CO; 

C + H 2 0 —> CO + H 2 ; 
CO + H 2 0 —» C0 2 + H 2 ; 

Tr 


t ▼ 

ASH CO, C0 2 , H 2 , H 2 0, CH 4 , C 2 H 2 , N 2 

Fig. 1 . Three zone equilibrium and kinetic free model of the downdraft gasifier. 


while Mott and Spooner [22] has assumed (2/3) O to be associated 
with hydrogen and (1/3) O with C. 

The association of fuel oxygen with carbon and hydrogen both 
seems to be logical if one studies both HHV predictions based on 
pure pertinent combustion reactions [14] and the devolatilization 
results [5,6,9,10,20,21,23]. In the present model of pyrolysis zone, 
based on these realizations and the fact that the affinity between H 
and O is much higher than that of C and O, it is assumed that 80% of 
fuel oxygen, i.e., (4/5) O is associated with fuel hydrogen in the form 
of H 2 O while, 20% of fuel oxygen, i.e., (1/5) O is associated with fuel 
carbon on mole basis and releases as CO and CO 2 on decomposition 
and the ratio of moles of CO and CO 2 is assumed to be inversely 
related with their molecular mass [5-8,13]. Fifty percent of 
balanced hydrogen in fuel releases as H 2 on decomposition [6,13] 
and the remaining 50% of balanced hydrogen in the fuel is released 
in the form of CH 4 and C 2 H 2 and the ratio of moles of CH 4 and C 2 H 2 
is assumed to be inversely related with their molecular mass 
[5-8,13]. The gas composition and temperature in this zone is 
governed by mass and energy balance. 

The product from the pyrolysis zone enters oxidation zone 
where incoming air reacts with these products and H 2 is fully 
oxidized to H 2 O [2,24-29]. Balanced oxygen will be consumed in 
char oxidation [2,3,27-31 ] due to a very high reactivity of char. The 
char oxidation yields CO and C0 2 , which are assumed to be 
inversely proportional to the exothermicity (nco/n<:o 2 = 3.5606) of 


The model formulated in this work consists of three constituent 
sub-models, one for each zone, i.e., pyrolysis, oxidation and 
reduction zone. The basic assumptions on which the overall model 
has been formulated are similar to two zone equilibrium and 
kinetic free model [ 1 ]. 

(1) The char is modeled as graphitic carbon [2,38,39]. 

(2) The process of gasification is isobaric [3,39,41]. 

(3) Char carry over from the gasifier is taken as 0% on molar basis 
[38,39]. 

(4) There is no heat transfer between the various zones [1,2]. 

(5) The overall direct heat loss from the gasifier is assumed to be 
related to both the equivalence ratio and the HHV of biomass 
material [2]. Heat loss depends on temperature levels of 
gasifier. The temperature levels in gasifier are governed by 
equivalence ratio and HHV of biomass. For gasification process 
requiring oxygen deficient atmosphere, the reactor tempera¬ 
ture increases with increasing equivalence ratio and HHV of 
feed material, therefore, the overall heat loss is assumed to be 
the product of equivalence ratio and HHV. For biomass having 
HHV in the range of 15-20 MJ/kg and operating in equivalence 
ratio range of 0.25-0.45, the overall heat loss will vary in the 
range of 3.0-6.0% [2,43]. Taking this into account the overall 
heat loss is assumed to be 10 % of the product of equivalence 
ratio and HHV [2,41-43]. This loss is distributed in pyrolysis, 
oxidation and reduction zones based on the study of temper¬ 
ature profiles of gasifier by Bliek [43] and Channiwala [2] 
where the maximum temperature is observed at oxidation 
zone, slightly lower levels in the reduction zone and still lower 
levels in pyrolysis zone. Typical average temperature is in the 
order of 850, 650 and 500 °C in oxidation, reduction and 
pyrolysis zone, respectively. Accordingly the heat losses in 
pyrolysis, oxidation and reduction zone are assumed to be 25, 
40 and 35% of this loss, respectively [2]. 

3.2. Formulation of pyrolysis zone model 
Assumptions: 

(1) 4/5 of the fuel oxygen is associated with fuel “h” in the form of 
H 2 0 [14-19,22]. 
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(2) 1/5 of the fuel oxygen is associated with fuel carbon and 
releases as CO and C0 2 on decomposition [14-19,22]. 

(3) The ratio of moles of CO and C0 2 is inversely related with their 
molecular mass, i.e., C0/C0 2 = 44/28 [5-8,13]. 

(4) 50% of available hydrogen in fuel releases as H 2 on decompo¬ 
sition [5,13]. 

(5) Balanced 50% of available hydrogen in fuel is released in the 
form of CH 4 and C 2 H 2 [5,13]. 

( 6 ) The ratio of moles of CH 4 and C 2 H 2 is inversely related to with 
their molecular mass, i.e., CH 4 /C 2 H 2 = 26.016/16.032 [5-8,13]. 

(7) The products of pyrolysis zone are C, CO, C0 2 , H 2 , H 2 0, CH 4 and 
C 2 H 2 [2-8]. 


Reaction stoichiometry: 
n c,p*C + n co 2 ,p*C0 2 + n CO p‘ CO + n C n 4 p«CH 4 

+ n H 2 ,p*H 2 + Hc 2 h 2 ,p*C 2 H 2 + n H 2 o,p*H 2 0 

+a(0 2 + 3.76N 2 ) —►Hqoxd’C + n co,oxD*CO (8) 

+ n co 2 ,oxD*C0 2 + ncH 4 ,oxD’CH 4 

+ n c 2 H 2 ,oxD , C 2 H 2 + Hh 2 o, oxd’H 2 0 + n N2; oxD*N 2 

Constituent balance: 

ttlCoXD = n C,P + n CO,P + n C0 2 ,P + n CH 4 ,P + 2«n C2 H 2 ,P 
= n C,OXD + n CO,OXD + n C0 2 ,OXD + n CH 4 ,OXD + 2n C2 H 2 , 0 XD ’ 


Reaction stoichiometry: 

CiH h 0 0 + wH 2 0 1 —►n C;P -C + nco 2 ,p’C0 2 + Hco,p*CO 
= ^ch 4 ,p*CH 4 + n H 2 ,p*H 2 + nc 2 H 2 ,p’C 2 H 2 + n H 2 o,p*H 2 0 (1) 
Constituent balance: 


ttlhoxD = 2 hh 2 ,p + 2hh 2 o,p + 4h C h 4 ,p + 2«n C2 H 2 ,p 
= 2-n C2 H 2 oxd + 4n CH4 oxd + 2n H20 ,oxD 

ttl°oxD = n co,p + 2n C o 25 p + h H 2 o,p + 2 a 
= n CO,OXD + 2n C 0 2? 0XD + n H 2 0,OXD 


( 10 ) 

( 11 ) 


ttlcp = 1 = ric p + n C o,p + n co 2 ,p + n CH 4 ,p + 2n C 2 H 2 ,p 
ttlhp = h + 2 w = 2 n H2 ,p + 2 n H 2 o,p + 4n CH4 ,p + 2 n C 2 H 2 ,p 

ttlop = o + w = n C o,p + 2 n C o 2 ,p + ^h 2 o,p 
Energy balance: 


( 2 ) 


(3) 


(4) 


tthioxD = (3.76*2)a = 2n N2 0 xd 
E nergy balance: 


E ni <T h°f + (^T,P - ^298) 


K,P 


= E n i^XD h°f + (hj OXD - ^298) 


j,OXD 




l,OXD 


( 12 ) 


(13) 


h°f + (h Tr - /i 298 ) n +w[/iQ f + (h Tr - h 2 9 g)]i o 


BM 


E n jT + (frr,P - ^ 298 ) . p+^l,P 



3.3. Formulation of reduction zone model 
Assumptions: 


3.2. Formulation of oxidation zone model 
Assumptions: 

(1) Hydrogen from the pyrolysis zone is fully oxidized to H 2 0 due 
to its highest burning velocity [2,24-29]. 

(2) The balance oxygen will be consumed in char oxidation due to 
much larger reaction area available for 0 2 adsorption on highly 
reactive pyrolysis char [2,3,27-31]. 

(3) CO, C0 2 distribution is considered to be inversely proportional 
to the exothermicity of their reactions, i.e., less exothermic 
reaction is assumed to be proceeding at a faster rate than 
higher one [2,3,27-29]. The principal char oxidation reactions 
considered are: 

C + lo 2 -CO (JHr = -110.6 kj/mole) ( 6 ) 

C + 0 2 ^C0 2 ( 2 lH R = —393.8 kj /mole) (7) 

These reactions will proceed in the inverse ratio of their 

exothermicity, i.e., nco,o/nco 2 ,o = 3.5606. 

(4) The CO and C0 2 from the pyrolysis zone is assumed to be 
simply added up to the CO and C0 2 produced by char oxidation. 

(5) CH 4 and C 2 H 2 are assumed to be carried forward to reduction 
zone due to their low burning velocities and lack of oxygen 
[2,27-29]. 

( 6 ) The products of the oxidation zone are char, CO, C0 2 , H 2 0, CH 4 , 
C 2 H 2 and N 2 [2-4]. 


(1) The products of the reduction zone are assumed to be CO, C0 2 , 
H 2 , H 2 0, CH 4 , C 2 H 2 and N 2 [2,5-8]. 

(2) CH 4 and C 2 H 2 from the oxidation zone are carried forward to 
product gas due to oxygen deficiency and consumption of char 
by reduction reactions [2,3,27-29,40]. 

(3) The gas phase homogeneous shift reaction is assumed to have 
a degree of reaction as a unity at the entry of reduction zone [ 2 - 
4,32-37]. 

C0 + H 2 0^C0 2 + H 2 (14) 

n CO,OXD • CO + H C o OXD • H 2 O -► n co OXD • C0 2 + H C o OXD • H 2 

(15) 

(4) The principal reduction reactions are Boudouard reaction and 
water-gas reaction. Both these heterogeneous reactions are 
assumed to consume char in the inverse ratio of their endo- 
thermicity [3,27-29]. This means that water-gas reaction 
proceeds faster than Boudouard reaction. A total of43.22% of char 
from the oxidation zone is consumed by Boudouard reaction and 
56.78% of char from the oxidation zone is consumed by water-gas 
reaction. Thus, both these reactions may be represented as: 

Boudouard reaction ; C + C0 2 ^2CO 
(JH r = 172.6 kj/mole) (16) 

0.4322n c ,oxD •mole + 0.4322n C oxD*CO 2 
—2*0.4322n CjO xD • CO 


( 17 ) 
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Fig. 2. Overall solution algorithm for a three zone equilibrium and kinetic free model [54]. 


Water-gas reaction; C + H 2 O^CO + H 2 (AH R 
= 131.4 kj/mole) (18) 


0.5678n c oxd’C H - 0.5678n^oxD*H 2 O—► 0.5678/1 coxd*^'O 
+ 0.5678n c oxd’H 2 (19) 

Reaction stoichiometry: 


n c,oxD*C + n co,oxD*CO + n C o 2 ,oxD , C0 2 + ncH 4 ,oxD*CH 4 
+ n H 2 o,oxD*H 2 0 + h C2H2 qxd.C 2H 2 + n Nz qxd*N 2 ->n c ^«C 

+ n co,R’CO + H C o 2 ,tf*C0 2 + n CH 4 ,R*CH 4 + n H2 ojr H 2 0 
+ n c 2 H 2 ,R* c 2 H 2 + n H 2 ,R‘H 2 + h N 2 r*N 2 
Constituent balance: 

ttlc^ = n c OXD + n CO,OXD + n C0 2 ,0XD + n CH 4 ,OXD + 2n C2 H 2j0 XD 
= n C,R + n CO,R + n C0 2 ,R + n CH 4 ,R + 2n C2 H 2 ,R 


( 20 ) 



[ R 


2n H; 


^h 2 ,oxd + 2hh 2 o,oxd + 4n C H 4 ,oxD + 2n C2 H 2 ,oxD 
2hh 2 ,r + 2hh 2 o,r + 4n C H 4? R + 2 /i C2 h 2 ,r 


Table 1 

Details of biomass with input parameters selected from published literature for 
model validation [4,50-53] 


Material/reference/molecular 

formula 

C/H molar 
ratio 

Equivalence 

ratio 

Moisture 
content(%) 

HHV 

(MJ/kg) 

Peach pits No. 2 [4] 

ClH 1 . 562 O 0.568 

0.640 

0.258 

10.9 

23.179 

Peach pits No. 3 [4] 

ClH 1 . 562 O 0.568 

0.640 

0.2638 

10.9 

23.179 

Peach pits No. 4 [4] 

ClH 1 . 562 O 0.568 

0.640 

0.2673 

10.9 

23.179 

Douglas fir wood blocks [4] 
C 1 H 1 . 434 O 0.594 

0.697 

0.2936 

5.4 

21.296 

Whole log wood 
chips No. 1 [4] 

C 1 H 1 . 4 g 25 O 0.713 

0.675 

0.2764 

18 

19.976 

Whole log wood chips No. 2 
[4] 

C 1 H 1 . 413 O 0.670 

0.708 

0.3322 

32 

19.976 

Whole log wood chips No. 5 
[4] 

C 1 H 1 . 413 O 0.670 

0.708 

0.2464 

9.80 

19.976 

Corn cobs 6" choke [4] 
ClH 1 . 541 O 0.782 

0.649 

0.2893 

11 

19.232 

Corn cobs 6" choke [4] 
ClH 1 . 541 O 0.782 

0.649 

0.2669 

11 

18.744 

Cubed alfalfa seed straw [4] 
C 1 H 1 . 572 O 0 . 78 i 

0.636 

0.2600 

7.90 

18.360 

Cubed alfalfa seed straw [4] 
C 1 H 1 . 572 O 0 . 78 i 

0.636 

0.245 

7.90 

18.360 

Birch wood [50] 

ClH 1 . 476 O 0.73 

0.678 

0.3268 

10 

18.900 

Coir dust [50] 

ClHi.28lOo.605 

0.781 

0.3135 

15 

19.800 

Subabul wood [51] 
ClH 1 . 451 O 0.697 

0.689 

0.317 

11 

20.017 

Cubed California 
cotton gin trash [52] 
ClH 1 . 582 O 0.689 

0.632 

0.3016 

23.50 

20.227 

Wood R. No. 98 [53] 

ClHi.4640o.695 

0.683 

0.323 

17.81 

19.668 

Wood R. No. 98 [53] 
ClH 1 . 464 O 0.695 

0.683 

0.277 

11 

19.668 

Wood chips R. No. 98 [53] 

ClH 1 . 464 O 0.695 

0.683 

0.2839 

16.52 

19.668 

Wood chips R. No. 920 [53] 
ClH 1 . 464 O 0.695 

0.683 

0.3251 

6.01 

19.668 

Wood chips R. No. 910 [53] 
ClHi.4640o.695 

0.683 

0.2778 

6.99 

19.668 

Wood chips R. No. 910 [53] 

ClH 1 . 464 O 0.695 

0.683 

0.285 

17.81 

19.668 

Wood chips R. No. 929 [53] 
ClH 1 . 464 O 0.695 

0.683 

0.2667 

7.41 

19.668 

Wood chips R. No. 1117 [53] 

ClH 1 . 464 O 0.695 

0.683 

0.2892 

11.26 

19.668 

Wood chips R. No. 910 [53] 

ClHi.4640o.695 

0.683 

0.2978 

16.52 

19.826 


ttloR = n co oxd + 2n C o 2? oxD + n H 2 o,oxD 
= n co,R + 2 h C o 2 ,r + Hh 2 o,r 

ttlnR = 2n N2 oxd = 2n N2 R 
Energy balance: 


Y) "K, 

k 


OXD 


h°f + (b T OXD - ^298) 


K,OXD 


E U i R + (^T,R - ^298) . r +^1,R 


3.4. Solution algorithm 


(23) 

(24) 



( 22 ) 


Independent solution of each zone of model is obtained by the 
method of iteration and coupled with each other to yield an overall 
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[Range of Input Parameters : ER = 0.245 to 0.3322, 

MC = 5.4% to 23.50%, HHV = 18.360 MJ/kg to 23.179 MJ/kg & 
C/H Molar Ratio = 0.632 to 0.781] 



- CO PREDICTED □ C02 PREDICTED 
A H2 PREDICTED O CH4 PREDICTED 
+ N2 PREDICTED 


AOWAE = 18.37% 
AOWBE = -0.55% 


Fig. 3. Model validation with published experimental results of the final gas composition. 


solution of the model. Fig. 2 depicts the overall solution algorithm 
of the model. The first zone is pyrolysis zone. Output of pyrolysis 
zone model becomes input for oxidation zone model and the 
output of oxidation zone model becomes input for reduction zone 
model. The computer program to solve each zone model has been 
developed in Turbo-C/C ++ language and the complete program 
interlinks all of the three zone solutions for getting the gas 
compositions, temperature at each zone and gasifier performance 
parameters. 

4. Model validation 

To establish utility of the present three zone equilibrium and 
kinetic free model for the complete range of solid feed stocks and 
operating conditions, the sufficient numbers (24 different biomass) 
of published data points have been selected for model validation. 
The various biomasses with input parameters selected are given in 
Table 1. In all cases char is modeled as graphitic carbon [2,38,39] 
balance char in reduction zone is assumed to be zero [38,39] and 
heat loss is assumed to be 10% of the product of equivalence ratio 
and HHV [2,41-43]. 

Figs. 3 and 4 shows the model validation with published 
experimental results of the final gas composition and LHV (lower 
heating value) of gas, respectively. C 2 H 2 is not compared due to 
absence of published experimental data. It can be seen from Figs. 3 
and 4 that the model predictions are relatively in good agreement 
with published experimental results as reflected from their error 
levels. The overall weighted absolute (AOWAE) and bias error 
(AOWBE) of the model over 24, published experimental data for gas 
composition is 18.37 and -0.55%, respectively, as shown in Fig. 3. 
The average absolute error (AAE) and average bias error (ABE) for 
LHV is 7.82 and 1.37%, respectively, as shown in Fig. 4. These results 
definitely, prove the worthiness of the model looking to the large 


numbers of variables involved in the published experimental data. 
It can therefore be stated that the present model developed with 
the simplest approach can be used to generate biomass specific 
performance data, which is needed for reactor design. 

5. Model predictions 

Main aim of any model lies in its utility as a key for design, scale- 
up and performance predictions. Accordingly, the present three 
zone equilibrium and kinetic free model, which has been validated 
with published experimental results to a sufficient degree of 
accuracy, has been used to study the influence of equivalence ratio 
and moisture content on the zone wise gas composition, temper¬ 
ature and gasifier performance parameters. The influence of 
equivalence ratio and moisture content on zone wise gas compo¬ 
sition, temperature and gasifier performance parameters have been 
studied for Subabul wood as a feed stock, the details of which is 
presented in Table 2. The salient features of these predictions and 
observations are presented in the following paragraphs. 

5.2. Influence of equivalence ratio (ER) 

Table 3 shows the gas composition of the pyrolysis zone. No 
effect of equivalence ratio has been observed on the gas composi¬ 
tion in pyrolysis zone, which is due to the fact that the present 
model is developed neither considering reaction rates nor consid¬ 
ering equilibrium. The pyrolysis zone product and temperature are 
obtained simply through mass and energy balance with few justi¬ 
fying assumptions. This fact indicates the limitation of this model 
because pyrolysis is governed by the rate of heating and tempera¬ 
ture levels [15-19]. 

Fig. 5 shows the influence of equivalence ratio on the gas 
composition of oxidation zone. Increasing N 2 with increasing 


CO 

E 


~o 

0) 

+■» 

o 

'S 

0) 


Q_ 



[Range of Input Parameters: ER = 0.245 to 0.3322, MC = 5.4% to 23.5%, 
HHV = 18.360 MJ/kg to 23.179 MJ/kg & C/H Molar Ratio = 0.632 to 0.781] 



LHVg Measured, MJ/Nm 3 


AAE = 7.82% 
ABE = 1.37% 


Fig. 4. Model validation with published experimental results of LHV (lower heating value) of gas. 
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Table 2 

Range of various parameters selected for model prediction 


Biomass Subabul 

HHV 

Range of moisture 

Range of 

Heat loss 

wood 

(MJ/kg) 

content(%) 

equivalence ratio 

(%) 

LlH 1 . 4513 O 0.6986 

19.777 

0-38 

0.25-0.45 

10 


equivalence ratio is quite obvious, while CH 4 , C 2 H 2 and H 2 0 are 
observed to be decreasing. This may be explained on the basis of 
the fact that CH 4 , C 2 H 2 and H 2 0 generations depends on elemental 
composition of biomass which is constant for given biomass 
material and hence the relative compositions will obviously 
decrease with increasing air intake and equivalence ratio. The CO 
and C0 2 levels are observed to be increasing with increasing 
equivalence ratio. This is obviously due to more and more oxidation 
of char with increasing equivalence ratio. Furthermore, the rate of 
increase in CO is observed to be higher than that of C0 2 . This may be 
attributed to the fact that under oxygen deficient atmosphere and 
the presence of highly reactive char, the partial oxidation of char 
will be dominant [30,31,44-46]. 

Fig. 6 depicts the influence of equivalence ratio on the gas 
composition of reduction zone. The decreasing CO and increasing 
C0 2 with increasing equivalence ratio is certainly due to oxidation 
of CO at higher equivalence ratio. The concentration of both these 
constituents is nearly equal at equivalence ratio of 0.3. This signifies 
that the gasifier must be operated below this equivalence ratio if 
one desires to achieve better gas quality. The hydrogen concen¬ 
tration reduces at a slower rate as compared to H 2 0 concentration 
with increasing equivalence ratio. This may be explained by the fact 
that both these constituents are involved in water-gas shift equi¬ 
librium, which consumes CO and H 2 0 and produces H 2 and C0 2 [2- 
4,33-35,39,47]. The CF1 4 is observed to be nearly constant due to 
model assumption. 

Fig. 7 shows the influence of equivalence ratio on the temper¬ 
ature at each zone. Temperature of pyrolysis zone is remaining 
constant with increasing equivalence ratio while temperature of 
oxidation and reduction zones are increasing with increasing 
equivalence ratio. The rate of increase of temperature as a function 
of equivalence ratio is higher in reduction zone as compared to 
oxidation zone, because the products of the oxidation zone enters 
at a high temperature in reduction zone plus the amount of char 
needed for endothermic reduction reactions is reduced, because at 
higher equivalence ratio most of the char is utilized in oxidation 
zone (as shown in Fig. 5). This is more clearly observed at equiva¬ 
lence ratio of 0.45 where both oxidation and reduction tempera¬ 
tures become nearly equal [2-4,34,35,39,48,49]. 

Figs. 8 and 9 presents the variation of gas yield and LHV of gas as 
a function of equivalence ratio. Obviously, gas yield in oxidation and 
reduction zones are increasing with increasing moles of air intake 


Table 3 

Gas composition of pyrolysis zone 


CO.p (%) 

co 2jP (%) 

H 2iP (%) 

H 2 O iP (%) 

CH 4 , p (%) 

C 2 H 2jP (%) 

C, P (%) 

T.p (K) 

5.94 

3.78 

8.05 

77.26 

3.08 

1.89 

17.18 

816.96 


while, gas LHV sharply decreases with increase in equivalence ratio, 
due to oxidation of gas phase constituents [2,4,34,35,39]. LHV of gas 
is optimum at equivalence ratio in the vicinity of 0.25 for woody 
biomass where solid carbon conversion is completed [2,3,4,34,41]. 
Fig. 10 shows the energy efficiency (which is the product of gas 
yield and gas LHV) of gasifier at reduction zone as a function of 
equivalence ratio where energy efficiency of gasifier is decreasing 
with increasing equivalence ratio, which is quite obvious. 

Thus, the overall effect of increasing equivalence ratio is to 
reduce the level of CO, H 2 and CH 4 and thereby to reduce the gas 
quality. These observations are quite consistent with the field 
observations as well as the predictions from the thermodynamic 
equilibrium and KF models [2-4,39,48,49] and thus once again 
validate the model qualitatively. 

5.2. Influence of moisture content (MC) 

Figs. 11-17 presents the influence of moisture content on the 
gas composition, temperature level and gasifier performance 
parameters at each zone with constant equivalence ratio of 0.3. 
The heat loss is assumed to be 10% of the product of HHV and 
equivalence ratio [2,41-43] and the char carry over is assumed to 
be 0% [38,39]. 

Figs. 11-13 represents the influence of moisture content on 
the gas composition at each zone. Fig. 11 shows the gas compo¬ 
sition of pyrolysis zone as a function of moisture content. It is 
observed that all concentrations except H 2 0 are decreasing with 
increasing moles of water vapor, which is quite obvious. Fig. 12 
shows the gas composition of oxidation zone as a function of 
moisture content. It is observed that all concentrations are 
decreasing with increasing moles of moisture except H 2 0. For the 
same biomass and constant equivalence ratio, these predictions 
are quite logical. Fig. 13 shows the influence of moisture content 
on final gas composition. Once again it is observed that CO, H 2 
and CH 4 concentrations are decreasing with increasing moisture 
content, which is obviously due to increased availability of moles 
of moisture in the system. 

Fig. 14 shows the influence of moisture content on the 
temperature at each zone. Decreasing trends of temperature at 
each zone as a function of moisture content is quite reasonable 
because higher moisture content consumes more amount of heat 
energy as a latent and sensible heat and thereby reduces the 
temperature levels [2,3,39,47]. 


[Material : Subabul Wood, MC = 15%, Heat Loss = 10%] 



■0— COoxd □ C02oxd X CH4oxd A C2H2oxd X H2Qoxd —©—N2oxd —I— Coxd 


Fig. 5. Influence of equivalence ratio (ER) on gas composition of the oxidation zone. 
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COr -eh- C02r -a- H2r H20r CH4r -©- C2H2r -h- N2r 


Fig. 6. Influence of equivalence ratio on gas composition of the reduction zone. 


[Material: Subabul Wood, MC = 15%, Heat Loss = 10%] 
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[Material: Subabul Wood, MC = 15%, Heat Loss = 10%] 
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Fig. 10. Influence of equivalence ratio on energy efficiency of the gasifier (reduction 
zone). 



Fig. 11. Influence of moisture content (MC) on gas composition of the pyrolysis zone. 


Fig. 7. Influence of equivalence ratio on temperature level at each zone. 



Fig. 8. Influence of equivalence ratio on gas yield of oxidation and reduction zones. 



[Material: Subabul Wood, ER = 0.3, Heat Loss = 10%] 



COoxd —b— C02oxd —H20oxd 
CH4oxd —l— C2H2oxd —©— N2oxd 
—i— Coxd 


Fig. 12. Influence of moisture content on gas composition of the oxidation zone. 

Figs. 15-17 represents the variations of gas yield, LHV of gas and 
energy efficiency (i.e., the product of gas yield and gas LHV) of the 
gasifier as a function of moisture content, respectively. Increasing 
gas yield and decreasing the trend of LHV with increase in moisture 
can easily be linked with the increased availability of moles of 
moisture in the system [2,4,39,47]. The overall energy efficiency at 
reduction zone does not show significant variation with increasing 
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[Material: Subabul Wood, ER = 0.3, Heat Loss = 10%] 



COr -B- C02r —A— H2r 
H20r CH4r —e— C2H2r 

—I— N2r 


Fig. 13. Influence of moisture content on gas composition of the reduction zone. 


[Material: Subabul Wood, ER = 0.3, Heat Loss = 10%] 
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[Material: Subabul Wood, ER = 0.3, Heat Loss = 10%] 
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Fig. 14. Influence of moisture content on the temperature level at each zone. 



0 — Gyieldp □ Gyieldoxd — A Gyieldr 


Fig. 15. Influence of moisture content on the gas yield at each zone. 

moisture content because it is simply the product of gas yield and 
gas LHV. Yield of gas increases while LHV of gas decreases and 
hence energy efficiency of gasifier does not show significant vari¬ 
ations. However, it must be noted that at very high levels of 


Fig. 17. Influence of moisture content on energy efficiency of the gasifier (reduction 
zone). 

moisture, reaction temperatures may not be favorable which will 
result in much inferior gas quality and hence efficiency. Thus it may 
be stated that the model very logically handles the predictions of 
gas composition and temperature levels in each zone as a function 
of moisture content, too. 

6. Conclusions 

Based on the formulation of the three zone equilibrium and 
kinetic free model of biomass gasifier, validation and the predic¬ 
tions of the following conclusions may be derived. 

(1) The comparison of the model predictions with published 
experimental results and illustrated closeness between the two 
within ±20% proves the quantitative and qualitative validity of 
the present three zone equilibrium and kinetic free model for 
most parameters of interest. 

(2) The present three zone equilibrium and kinetic free model is 
the simplest and offers a novel approach of modeling the 
gasifier. This model clearly provides the operating range of 
equivalence ratio and moisture content for the woody biomass 
materials. Furthermore, this model facilitates the prediction of 
maximum temperature in oxidation zone of gasifier, which 
provides useful information for the design of gasifier and the 
selection of material for the construction. 

(3) The merits of the model lies in the fact that it is capable of 
handling predictions for all category of biomass materials with 
wide operating range of equivalence ratio and moisture 
content in all the three principal zones of gasifier. 
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